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We theoretically examine Raman photoassociation of a Bose-Einstein condensate, revisiting stim- 
ulated Raman adiabatic passage (STIRAP). Due to collisional mean-field shifts, efficient molecular 
conversion requires strong coupling and low density, either of which can bring about rogue photodis- 
sociation to noncondensate modes. We demonstrate explicitly that rogue transitions are negligible 
for low excited-state fractions and photodissociation that is slower than the STIRAP timescale. 
Moreover, we derive a reduced-parameter model of collisions, and thereby find that a gain in the 
molecular conversion efficiency can be obtained by adjusting the atom-atom scattering length with 
off-resonant magnetoassociation. This gain saturates when the atom-atom scattering length is tuned 
to a specific fraction of either the molecule-molecule or atom-molecule scattering length. We con- 
clude that a fully-optimized STIRAP scheme may offer the best chance for achieving coherent 
conversion from an atomic to a molecular condensate with photoassociation. 

PACS numbers: 03.75.Fi,05.30.Fk,32.80.Wr 
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I. INTRODUCTION 

Photoassociation occurs when a pair of colliding ul- 
tracold atoms absorb a laser photon 0, thereby jump- 
ing from the two-atom continuum to a bound molecular 
state. If the initial atoms are Bose-condensed H B 0, H[ , 
then the subsequent molecules will also form a Bose- 
Einstein condensate H H 0- On the other hand, the 
process known as the Feshbach resonance 9], which we 
refer to as magnetoassociation, occurs when two ultra- 
cold atoms collide in the presence of a magnetic field, 
whereby the spin of one of the colliding atoms flips, and 
the pair might then jump from the two-atom continuum 
to a bound molecular state 

[liJIllI]. As with collective 
photoassociation, the so-formed molecules will comprise 
a Bose-Einstein condensate (BEC) if the incident atoms 
are themselves Bose-condensed 0, 0, Q . 

Early theories of collective association accounted only 
for the condensates, ne glecting noncondensate modes M, 
US III 13 Rogue THEi, or unwanted [HEEl, 
transitions to noncondensate modes occur because the 
dissociation of a zero-momentum BEC molecule need 
not take the atoms back to the zero-momentum atomic 
condensate, but may just as well end up creating two 
noncondensate atoms with equal-and-opposite momenta. 
Since the collective condensate coupling scales like the 
square root of the laser intensity (Feshbach-resonancc 
width) and the dissociation rate scales like the intensity 
(width), rogue dissociation is expected to play a dom- 
inant role in strong photoassociation (magnetoassocia- 
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tion) [H|T1|T1|T1 E3. 

Pioneering experiments [2I] with photoassociation 
of 87 Rb condensate were found to be just on the 
verge of coherent atom-molecule conversion. Next- 
generation Na [2l| and 7 Li [2I] experiments were aimed 
at the strongly interacting regime, and probed the pre- 
dicted 0, Qi| El photoassociation rate limit. Mean- 
while, groundbreaking experiments in magnetoassocia- 
tion of a Na condensate demonstrated a tunable scat- 
tering length [23, as well strong condensate losses for 
sweeps of the magnetic field across resonance [24| . Sub- 
sequent experiments in Feshbach-tuning the scattering 
length led to the formation of stable 85 Rb [2{| and 
137 Cs ;2f| condensates-which are otherwise incondensi- 
ble in significant fractions due to uncooperative (zero- 
field) scattering lengths (2^, as well as a controlled col- 
lapse of a condensate with bursts of atoms emanating 
from a remnant condensate 28], a counterintuitive de- 
crease in condensate losses for an increasing interaction 
time p9j . and collective burst-remnant oscillations |3C| . 
It turns out that rogue dissociation is at the heart of 
these magnetoassociation experiments 0, Isij . 
Most recently, short-lived quantum degenerate molecules 
have been created by sweeping a BEC across a Feshbach 
resonance |35| . 

When it comes down to creating a molecular conden- 
sate, the catch to photoassociation is that it generally 
occurs to an electronically-excited state, and the subse- 
quent irreversible losses to spontaneous decay defeat the 
purpose of molecular coherence. Adding a second laser 
to drive molecular population to a stable electronic state, 
stimulated Raman adiabatic passage (STIRAP) j^l| in 
photoassociation I37I l38l l39| of a Bose-Einstein conden- 
sate 0, ^1 ^1 1 1 jl | has been proposed as a means for 
avoiding radiative decay. Likewise, combining STIRAP 
with near-resonant magnetoassociation is also a viable 
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means to avoid spontaneous decay when creating a sta- 
ble molecular BEC @ |H |H although the laser 
spectroscopy (bound-bound-bound) is different from the 
usual photoassociation (free-bound-bound). 

The hallmark of STIRAP is the counterintuitive pulse 
sequence, which in photoassociation of a BEC amounts to 
adjusting the two lasers so that in the beginning, when 
the system is mainly atomic condensate, the strongest 
coupling is between the excited-molecular and stable- 
molecular condensates, while in the end, when effectively 
everything is in the target state, the strongest coupling 
is between the atomic and excited-molecular condensate. 
As the population is transferred between the atomic and 
stable-molecular condensates, the state with the larger 
population is always weakly coupled to the electronically- 
excited molecular condensate, and the subsequently low 
(ideally zero) population reduces (eliminates) radiative 
losses. In principle, STIRAP allows for complete conver- 
sion from an atomic to a molecular condensate. 

In practice, however, free-bound-bound STIRAP re- 
lies on a superposition that includes the atoms and 
stable molecules but excludes the electronically-excited 
molecules |40). and the conversion efficiency is reduced 
from unity when this so-called dark state 0] is disrupted 
by collisions between particles. Specifically, mean- 
field shifts due to collisions between condensate parti- 
cles (atom-atom, atom-molecule, and molecule-molecule) 
make it difficult to achieve STIRAP by moving the sys- 
tem off the required two-photon resonance [ill |42T | . Such 
two-photon frequency shifts lead to a larger excited-state 
fraction, and thus more losses to radiative decay, which 
of course reduces the molecular conversion efficiency, de- 
feating the purpose of collective conversion. Neverthe- 
less, the collisional interaction strength is directly propor- 
tional to density, and the conversion efficiency can be im- 
proved by treating the density as an optimization param- 
eter (Z^. Another collision-avoidance option is STIRAP 
in an optical lattice with two particles per site |49(, al- 
though this scheme is not many-body coherent. 

The purpose of this Article is to re-investigate Raman 
photoassociation of a nonideal Bose-Einstein condensate, 
focusing on improving the efficiency of STIRAP from an 
atomic to a molecular condensate. In Sec. [H] we intro- 
duce the model and briefly review the idea of improving 
the molecular conversion efficiency by reducing the den- 
sity. Because this system is in fact strongly coupled, for 
low density or otherwise, and because rogue photodis- 
sociation to noncondensate modes can be prominent for 
strongly coupled systems, Sec. Illll explicitlv demonstrates 
that rogue transitions are negligible when the excited- 
state fraction is low and photodissociation is slow com- 
pared to the STIRAP timescale. The main contribution 
of this work is given in Sec. IIVI which illustrates that 
the conversion efficiency can be further improved by tun- 
ing the atomic scattering length with off-resonant mag- 
netoassociation. Overall discussion is given in Sec. 
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FIG. 1: Three-level illustration of coherent free-bound- 
bound photoassociation, where N atoms have assumedly 
Bose-condensed into state The free-bound and bound- 
bound Rabi frequencies are k and f2, respectively. Similarly, 
the two-photon and intermediate detunings are A and 5. The 
wavy line denotes the irreversible losses that STIRAP is in- 
tended to manage. 



II. 



DENSITY TUNING STIRAP 



Turning to the situation of Fig. ^ we assume that N 
atoms have Bose-condensed into the same one-particle 
state, e.g., a plane wave with wave vector k = 0. Pho- 
toassociation then removes two atoms from this state 
|1), creating a molecule in the excited state |2). Includ- 
ing a second laser, bound-bound transitions remove ex- 
cited molecules from state \2) and create stable molecules 
in state |3). In second quantized notation, boson anni- 
hilation operators for atoms, primarily photoassociated 
molecules, and stable molecules are denoted, respectively, 
by a, b, and g. The laser-matter interactions that drive 
the atom-molecule and molecule-molecule transitions are 
characterized by their respective Rabi frequencies k and 
n. The two-photon and intermediate detunings are de- 
noted by A and S , and the spontaneous decay rate of the 
electronically-excited molecules is 7 S . 

We also introduce atom-atom, atom-molecule, and 
molecule-molecule contact interactions with respective 
coupling strengths: 
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(lb) 
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where ciy is the particle-particle scattering length, V is 
the quantization volume that is expedient in a particular 
context (e.g., cubic box or spherical cavity), and m is 
the mass of an atom. For the time being, we assume 
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that the spontaneous decay rate of the excited- molecular 
condensate is sufficiently large to justify neglect of the 
mean-field shifts for the electronically excited molecular 
state 1 2). For future reference, we define Ay = VAy. 

The three-mode Hamiltonian for this freely-interacting 
system can be written 



H 

7T 



= iAaW (5- \i ls )b% 

— \ [n^aa + a* cr 6) + fi(. 9 t 6 + 6 t <?)] 
+i(Ana t a t aa + X^g^gg) + X l3 a)g^ag. (2) 



The mean-field equations for collective two-color pho- 
toassociation of a freely-interacting gas are obtained from 
the Heisenberg equations of motion: 



m 
ib 

ig 



AnH 2 
b- 
A33I3 



(6 - ±i 7i 

(A 13 |a' 2 ' * '~' 2 



+ A 13 \g\ 2 )a- X a*b, 
\ (xaa + ng) , 



g - h nb - 



(3a) 
(3b) 
(3c) 



Here the effects of Bose enhancement |40| have been 
included by scaling the mean-field amplitudes accord- 
ing to x — > x/yN (where x = a,b,g). The mean- 
field probability \a\ 2 (|fr| 2 ,|g| 2 ) is then of the order 
of unity (half), and the coupling strengths have been 



redefined as \ 
The transient 

x(t) 



Nk and A, 



NX, 



■■ pXij. 

STIRAP pulses are taken as Gaus- 
= xooxp[-(t- D x f /T 2 } and tt(t) = 



sian, i.e. 

fi exp[- {t~D 2 ) 2 '/T 2 ]. 

Assuming an ideal gas (Ay = 0), it is easy to show 
from the steady-state limit {x = 0) that Eqs. © reduce 
to an algebraic system that can be plugged trivially into 
Mathematica yielding exact solutions 0- Here we 
scale the atomic (molecular) amplitude(s) asa-> e lAt ^ 2 a 



(b 



JAti 



JAt 



g), and keep only normalized solu- 

2\ 



9 ~* e , 

tions [\a\ z +2(|&| 2 + \g\ 2 ) = 1] that are real for real values 
of f2/x — > (0, 00). The dark state is given as 



A 
«o 
b 
,9o 



\ 




= 0, 




(4a) 
(4b) 
(4c) 
(4d) 



It is evident from Eqs. I0J that, for a two-photon res- 
onance (A = 0), the dark state is an eigenstate of the 
system with all atoms (stable molecules) for f2/x — * 
(fi/x — > 00), and zero excited molecules; hence, for an 
initial BEC and adiabatic counterintuitive pulses, the 
two-photon-resonant system will follow this eigenstate as 
it evolves from atoms into stable molecules. The presence 
of collisions, which introduce a time-dependent shift of 
the two-photon detuning, means that the system cannot 



start out in, nor subsequently follow, the dark state, and 
a nonzero excited-state fraction will ensue. 

To illustrate numerically, we concentrate on explicit 
parameters for a dilute quantum-degenerate ga s of 8 7 Rb 
atoms 42j: 7s = 7.4 X H^B" 1 , Xo = 2.1 x 10 6 Vp7/^s-\ 
po = 4.3 x 10 14 cm" 3 , An = 4.96 x Kr n cm 3 /s, A13 = 
—6.44 x 10~ n cm 3 /s; although unknown, the stable- 
molecule mean-field shift A33 — 2.48 x 10~ n cm 3 /s is 
estimated by assuming equal atom-atom and molecule- 
molecule scattering lengths. The effect of collisions is 
demonstrated explicitly in Fig. El Setting the mean-field- 
shift terms to zero, Ay = 0, we see in Fig. [21 (a) that 
STIRAP works exactly as expected for a near-adiabatic 
pulse sequence, losing only a small fraction of particles 
to spontaneous decay. Reinstating collisions, Fig. |2 (b) 
shows an order of magnitude increase in the excited-state 
fraction, and a corresponding decrease in the conversion 
efficiency of some two orders of magnitude, as sponta- 
neous decay gets the upper hand. 

Nevertheless, since the collisional mean-field shifts are 
directly proportional to the density (Ay tx p), then we 
expect that reducing the density of the initial BEC could 
improve the conversion efficiency of STIRAP from an 
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FIG. 2: Stimulated Raman adiabatic passage in photoasso- 
ciation of a nonideal Bose-Einstein condensate. For a given 
Xo(p), the pulse parameters are f2o = Xo, T = 5 x 10 3 /xo, 
Di = 4.5T and D 2 = 2.5T. The two-photon (intermediate) 
detuning is A = (S = xo)- Note that \g\ 2 = 1/2 is actually 
complete conversion, since two atoms are destroyed to cre- 
ate a molecule, (a) For an ideal gas with p — po, SFIRAP 
mostly defeats spontaneous decay, as expected for a moder- 
ately adiabatic pulse sequence, (b) Adding collisional inter- 
actions shifts the system off resonance, leading to a larger 
excited-state fraction and, thus, more losses to spontaneous 
decay, (c) Density-tuning the STIRAP efficiency: Losses to 
spontaneous decay decrease as collisional mean-field shifts are 
marginalized for lower densities, but the STIRAP timescale 
(oc 1/ yfp) increases as well, and spontaneous decay eventually 
regains its dominance. Here a larger peak bound-bound pulse 
was used, fio = 50xo, as discussed in the text. 
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atomic to a molecular condensate |43j. This idea was 
spawned in our development of Feshbach-stimulated pho- 
toproduction of stable molecular condensates (see 
also related work |5l|]). Consider for example the den- 
sity p = 4.3 x 10 12 cm -3 , so that xo = 2.1 x 10 5 s _1 , 
An = 213 s- 1 , A 33 = 107 s" 1 , and Ai 3 = 277 s" 1 . The 
mean-field shifts are then roughly three orders of mag- 
nitude smaller than the peak Bose-enhanced free-bound 
coupling xo, and since xo sets the timescale for atom- 
molecule STIRAP [4(J], we expect a smaller role for par- 
ticle interactions compared to when p = p . This intu- 
ition is confirmed in Fig. (c) . Note that pulses with 
asymmetric heights |4l|, |42( improve the short-pulse effi- 
ciency of STIRAP by more than an order of magnitude 
over pulses with symmetric heights, as can be seen by 
comparing the final stable fraction in Fig. [21 (b) with the 
stable fraction for p/p = 1 in Fig. |2] (c). The opti- 
mal density for these pulse parameters is p ~ 10 12 cm -3 , 
for which conversion to a molecular BEC occurs on a 
timescale T = 5 x 10 3 /xo ~ 50 ms. Optimizing the de- 
tunings and pulse parameters |42l s hould give even fur- 
ther improvement (see also Ref. |45jV 



III. ROGUE PHOTODISSOCIATION TO 
NONCONDENSATE MODES 

As mentioned in the introduction, transitions to non- 
condensate modes occur because, although formed from 
a pair of zero-momentum BEC atoms, a condensate 
molecule need not photodissociate back to the atomic 
condensate, but may just as well form a noncondensate 
atom pair with equal and opposite momentum. Although 
experiments remain inconclusive on this score |2lL 12^ , a 
rate limit is to be expected when the conversion to rogue- 
dissociated atom pairs dominates over the formation of 
molecular condensate |0, H3, [Til . Similarly, the losses 
for across- resonance sweeps |25|. as well as the counter- 
intuitive losses 29], observed in 85 Rb magnetoassociation 
experiments can be viewed as collective rapid adiabatic 
passage from BEC to rogue-dissociated atoms pairs (3^ ; 
also, the remnant-burst oscillations [3(j are identified as 
Ramsey fringes in the evolution between an atomic con- 
densate and a molecular condensate dressed by noncon- 
densate atom pairs 152, H3, 0] • 

The prime indicator of rogue dominance is the density- 
dependent frequency uj p = hp 2 / 3 /m 0,0,0. When 
the collective-enhanced atom-molecule coupling satisfies 



realistic model 16] of rogue dissociation to nonconden- 
sate modes in order to determine its impact on photoas- 
sociative STIRAP in a freely-interacting Bose-Einstein 
condensate. Accordingly, the Hamiltonian (J2J is adapted 
to read 

— = %Aata + 5tfb- | [{ntfaa + ng^b) +H.c] 
+\{\i 1 a i a)aa + \ 33 g> g^ gg) + \ 13 <Jg*ag 
+3 ^2 [( A + £k ) a k flk _ ( &tfl k a -k + H.c. 



(5) 



Here ek = hk 2 /m is the energy of a noncondensate pair, 
and /k describes the wavenumber (energy) dependence 
of the noncondensate coupling. In order to isolate the 
effects of rogue dissociation, we have dropped the spon- 
taneous decay term ij s b^b. Also, anticipating low popu- 
lations, we have neglected collisions with noncondensate 
atoms. The corresponding mean-field theory is derived 
from the Heisenberg equations of motion: 



ib = 



iC(e) 
iP(e) 



5b -\ 



[iA + Aii|a| 2 +Ai 3 M 2 ] a~ X a*b, (6a) 

Xa 2 + ng + dJde^f(e)C(e) , (6b) 

[Ai 3 |a| 2 + A 33 | 3 | 2 ] 3 -ir!6, (6c) 
[A + e]C(e)-^ef(e)[l + 2P(e)]b, (6d) 

f(e)[b*C{e)-C*(e)b]. (6e) 



(2™ p ) 3 / 2 



In Eqs. JEJ, we have taken the continuum limit 

1 



k 



A 2 3/2 



deG(e), 



(7) 



and we also have introduced the rogue and normal den- 
sities, C(e) = </£ (a k a_ k ) /(27njp /4 ) and P(e) = (a^Ok), 

as well as the rogue coupling £(t) = x(i)/(27njp^ 4 ). 

Before moving on, we discuss basic renormalization. 
For xo = = and ignoring the normal density, sim- 
ple Fourier analysis of an initial excited-bound molecular 
condensate gives the below-threshold binding energy as 
the real and negative pole of 



lu — S — £(o->) + it] = 0, 



(8) 



A'o 



> 



Lu p , transitions 



to 



, . , where ri = + and the molecular self-energy is defined as 

noncondensate modes are ex- ' OJ 



pected to dominate atom-molecule conversion. For the 
previous example of p = 4.3 x 10 12 cm -3 , we find xo ~ 
100 Up] moreover, even for a typical density (p = po = 
4.3 x 10 14 cm- 3 ), we find xo ~ 50a; p . It is therefore 
reasonable to suspect that rogue photodissociation could 
play a dominant role in STIRAP from an atomic to a 
molecular condensate, for low-density systems or other- 
wise. 

The purpose of this section is to apply our minimal- yet- 



E(w) 



l£2 
2 t - 



/ def{e) 



(9) 



The simplest energy dependence for the continuum is one 
that obeys the Wigner threshold law up to some abrupt 
cutoff: / 2 (e) = 0(cm — e )- The detuning (binding en- 
ergy) then picks up a term 2(0) = £ 2 J~^m ■ I n principle, 
the cutoff is infinite, and therefore so is the continuum 
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shift of the molecular binding energy. To account for 
this divergence, one defines the so-called physical detun- 
ing 5 — S — £ 2 y / £M, which is finite by definition in the 
limit of an infinite cutoff. 

In practice, any numerical procedure employs a finite 
cutoff, and the finite shift is accounted for in exactly 
the same manner. Hence, the intermediate detunings in 
Eqs. are taken as physical (renormalized) detunings. 
Anticipating low noncondensate fractions, the normal 
density-Bose enhancement of the noncondensate modes- 
was neglected. Reminiscent of our original nondegenerate 
quasicontinuum model of photoassociation [3^, we in- 
troduce N qc = 2 x f 3 quantum-degenerate quasicontin- 
uum states and a cutoff ejvf = 10 xo, sufficient to deliver 
convergence and keep numerical artifacts to a minimum, 
where convergence is determined by the requirement 

H 2 + 2[|6| 2 + |<?| 2 ]+ J de\C(e)\ 2 = l. (10) 

That said, we report numerical STIRAP solutions of 
the mean-field theory 10 for rogue transitions to non- 
condensate modes in Raman free-bound-bound photoas- 
sociation of a BEC. The focus is on the lowest density 
from Fig. |3 p = 4.3 x 10 10 cm~ 3 , as an explicit example 
in Fig. |3 and the key results are summarized in Tabled 
In this case, not only is the coupling strong to begin 
with, but we have an especially low density. Here the 
rogue coupling is a bit stronger, due to the low density, 
leading to a slightly larger noncondensate fraction. But, 
due to the fact that mean-field shifts are effectively elim- 
inated for the lowest density, the bound state fraction 
is in this case the lowest. In the end, even for this low 
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FIG. 3: Rogue photodissociation to noncondensate modes in 
stimulated Raman adiabatic passage in photoassociation of 
an interacting Bose-Einstein condensate. The STIRAP pa- 
rameters are 8 = S7 = Xo, T = 5 x 10 3 /xo, D\ = 4.5T, and 
D2 = 2.5T. The fraction of rogue (noncondensate) particles, 
f r = Jde |C(e)| 2 , is clearly negligible, and STIRAP runs even 
better than in the ideal case [see Fig|5|(a)] because photodis- 
sociation is slower than spontaneous decay. 



TABLE I: Overall role of noncondensate modes for stimu- 
lated Raman adiabatic passage in photoassociation of an in- 
teracting Bose-Einstein condensate. The STIRAP parame- 
ters are as in Figure |3 Recall that the density scale is po = 
4.3 x 10 14 cm -3 , the rogue indicator is uj p = hp 2 ^/m, and 
the fraction of noncondensate particles is f r = f de|C(e)| 2 . 
Whereas the system is strongly coupled even for large densi- 
ties, i.e., xo > <+>p, the excited-state fractions are low enough, 
and rogue dissociation slow enough, to guarantee that a neg- 
ligible fraction of the population is transferred to the noncon- 
densate modes. 
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0.49983 
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2.34 x 10 2 


5.43 x 10~ 7 


1.54 x lQ- 


10 


0.49952 



density, we see that rogue photodissociation to noncon- 
densate modes does not affect the final stable-molecule 
conversion efficiency. We have confirmed that, for arti- 
ficially large rogue couplings £ > 1, i.e., faster photodis- 
sociation, rogue noncondensate modes indeed begin to 
play a more invigorated role. Also, if the timescale for 
STIRAP is lengthened too much, we expect that rogies 
should likewise begin to dominate. 



IV. FESHBACH TUNING STIRAP 

As we have already seen, reducing the density of the 
initial atomic condensate acts to marginalize collisional 
mean-field shifts of the two-photon resonance, leading to 
an improvement in the efficiency of STIRAP from an 
atomic to a molecular condensate. With this idea in 
mind, it seems possible that collisional mean-field shifts 
could be further marginalized by adjusting the atom- 
atom scattering length via off-resonant magnetoassoci- 
ation. We address this problem by first rewriting the 
Hamiltonian J5J in a form that makes the mechanics of 
the Feshbach-tuned efficiency more transparent. 

Neglecting rogue photodissociation to noncondensate 
modes as per the previous section, consider now the total 
number of particles, which is given by the operator 

N = a^a + 2(tfb + g^g). (11) 

Due to the damping term (oc ij s ) in the Hamiltonian J2Jl, 
it cannot be said that N is a conserved quantity. Never- 
theless, it is certainly true that N commutes with the 
Hamiltonian, i.e., [H,N] = 0, regardless of the non- 
Hermitian nature of H. It is then reasonable to pre- 
sume that one may add multiples of N to the Hamilto- 
nian without altering the essential physics. For example, 
it is implicitly accepted that the two free Hamiltonians 
Ho/h = iAata + Stfb and H /h = -Agig +{6- A)tfb 
will give the same physics, and the two are of course re- 
lated by the addition (subtraction) of the term ^HAN 
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[see also derivation of Eqs. Q]. Similarly, it is also rea- 
sonable to expect that subtracting the term ^HXn N 2 will 
leave the dynamics unchanged, an operation that results 
in the Hamiltonian 

-\ [K(tfaa + aVo) + Q(g^b + tfg)] , (12) 

where the parameters A 33 = A33 — 4An and A' 13 = 
A13 — 2An now characterize the mean-field shifts due to 
collisional interactions, and we continue to neglect the 
mean-field shifts associated with electronically-excited 
molecules. 

The reduced-parameter Hamiltonian (|12J) gives the 
mean-field equations 

1a = (±A + A' 13 \g\ 2 )a- X a*b, (13a) 
ib = (6 - |t7») b - \ { X aa + Qg) , (13b) 
ig = (A' 13 | a | 2 +A3 3 | 5 | 2 ) ff -i™. (13c) 

We have made an extensive numerical comparison be- 
tween the solutions to reduced-parameter mean-field 
model (|13|l and the solutions to the conventional mean- 
field model (JSJ). As expected, the conventional and 
reduced-parameter mean-field equations of motion lead 
to exactly the same physics. Moreover, we have also con- 
firmed that, barring a fluke in the scattering length(s), 
collisions with the excited-state are safely ignored. 

Besides reducing the number of parameters, it seems 
that A 33 = 0, which corresponds to an = 033/8, is a spe- 
cial case that could further improve the molecular con- 
version efficiency. Indeed, reducing the scattering length 
produces another marked gain in the molecular conver- 
sion efficiency [Fig. 01 (a)] . The gain saturates because 
An is steadily decreasing, so that A 33 — » A33. Similar 
saturation is shown in Fig. 0] (b) for A'i 3 = 0, but any 
gain here is outweighed by the fact that an < desta- 
bilizes the initial BEC against collapse, although a sys- 
tem where 013 > may prove otherwise. Again, further 
improvements can be achieved by optimizing the laser 
parameters 0>EH- 

Before closing, we extend our model to explicitly de- 
scribe tuning the scattering length via magnetoassocia- 
tion. Hence, in addition to the primary-photoassociation 
molecular state, the BEC atoms are coupled to a sec- 
ond molecular state, denoted by the operator 0, with a 
Feshbach resonance. For well-resolved resonances, light- 
driven transitions involving the Feshbach state are ne- 
glected just as, for example, direct free-bound photoasso- 
ciation transitions are neglected in Feshbach-stimulated 
photoproduction of molecular BEC |4J, |4(j . The Hamil- 
tonian is then appended to read 

H a 

— = ...+ Q B <P<j) = (c/> 1 aa + a 1 a t <j>) , (14) 

n v N 

where w B = [uj b — i(T d + r p a'a)]. Here T d is the rate 
of magnetodissociation to noncondensate modes and T p 



c 0.35 
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FIG. 4: Feshbach-tuning the efficiency of stimulated Raman 
adiabatic passage in photoassociation of an interacting con- 
densate. Here the density is p — 4.3 x 10 12 cm -3 , and the 
STIRAP parameters are the same as in Fig. |5] (c) for asym- 
metric pulses, (a) Compared to native collisions (an = 033), 
a Feshbach-tuned atomic scattering length offers additional 
efficiency gain, which saturates at an = 033/8. (b) Similar 
gain saturation also occurs at an = 3ai3/8 < 0. 

is the loss rate due to collision-induced vibrational re- 
laxation; also, the magnetoassociation coupling is a = 
y / 2i:p\aii\A p A B /m, and the energy (detuning) of the 
magnetoassociation state is hw B = sgn[an]A„(S — B a ), 
with A^ the difference between the atom pair and 
molecule magnetic moments and is the width of the 
resonance. Anticipating lowest-order perturbation the- 
ory, we have neglected molecule-molecule collisions. 

When the system is off-resonance (io B 3> a,Td, pT p ), 
it is easily shown from the Heisenberg equation of mo- 
tion for the molecular operator that, to lowest order, off- 
resonant magnetoassociation can be included simply by 
making the substitution 



An = An 


1 - 


,(T d 
-tl 










LO B J _ 



According to Fig. we are looking at about an order 
of magnitude decrease in the scattering length. For con- 
creteness, we assume SJ[An] = An/10, which translates 
into the magnetic field position B = Bq + 0.9 A B . The 
question is then whether this is sufficiently detuned from 
resonance to justify neglect of dissociative and relax- 
ational losses. The ultracold magnetodissociation rate 
is 52] Td = auA p A B p/h 2 , where p is the relative atom- 
pair momentum; hence, using the uncertainty relation 
p-V*p ~ h/2, we find T d /uj B ~ ail £^A B /(B - B ) ~ 
1(T 3 . Estimating T p ~ 10~ 10 cm 3 /s H2, A B ~ 0.1 G and 
A p = p B (i.e., a Bohr magneton), then the loss term due 
to vibrational quenching is roughly pT p /ui B ~ 10~ 4 for 
p ~ 10 12 cm~ 3 . It should therefore be possible, using 
one of the more than forty Feshbach resonances found in 
87 Rb 53], to improve the STIRAP efficiency by tuning 
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the atom-atom scattering length. 

V. DISCUSSION 

Admittedly, we have neglected an explicit trapping po- 
tential for either the atoms or the molecules, a move 
which we now justify. From Fig. (c), the timescale for 
STIRAP is T = 5x 10 3 /xo ~ 24 ms. For N = 5x 10 5 con- 
densate atoms in a spherically symmetric trap, let us say 
that the reduction in density from po = 4.3 x 10 14 cm~ 3 
to p ~ 10 12 cm -3 is achieved solely by reducing the fre- 
quency of the trap from v t — 100 s^ 1 to Vt ~ 1 s _1 . The 
STIRAP timescale is then fast enough to allow inclusion 
of the trap with an average over a Thomas-Fermi den- 
sity distribution which should not drastically alter 
the predicted optimal density. 

Also, the Raman-formed molecules can be vibra- 
tionally hot |42j . and molecule-molecule and atom- 
molecule collisions may thus foster relaxation to lower- 
lying vibrational levels. This vibrational quenching is 
accounted for in exactly the same fashion as already pre- 
sented in Sec.0]for the Feshbach molecular state. An ex- 
act value for the 87 Rb2 quenching rate (per unit density) 
is unknown (hence the neglect in Ref. |43), although the 
upper bound T' p < 10~ 10 cm 3 /s has been measured |20| . 
At this rate, vibrational quenching could well be an issue, 
even for the optimal density p ~ 10 12 cm~ 3 . However, 
a lower peak bound-bound Rabi coupling (CIq/xq = 50 
compared to 10 3 [42]) means that, for the same bound- 
bound laser intensities, it should be possible to target 
lower-lying levels. Because the vibrational-relaxation 
rate decreases with decreasing binding energy |54j |. tar- 
geting lower-lying vibrational levels should also lead to 
slower relaxation rates. In the end, vibrational quenching 
should be manageable. 

The remaining concern is whether the photoassociation 
laser will overlap with the expanded BEC cloud for user- 
friendly intensities. If not, then the cloudsize could be 
reduced by adjusting the number of condensate atoms; 
additionally, Feshbach-tuning the scattering length will 
also reduce the cloudsize. In particular, if the density 
p = 4.6 x 10 12 cm -3 p is achieved by lowering the 
number of atoms to N = 10 4 , the scattering length to 
&il = an/8, and the trap frequency to v t — 3s _1 , then 
we find that the cloudsize is only barely increased from 



Rq = 8.9 /im to R — 11 /im. By optimizing the number of 
atoms, the trap frequency and atomic scattering length, 
overlap between the BEC cloud and the photoassociation 
laser should be manageable for user-friendly intensities. 

Previously [15J , we have argued that experimental Ra- 
man photoassociation of a BEC (2(j is on the verge of 
coherent conversion. The experiments in question are 
performed with continuous-wave lasers and spontaneous 
decay from the primary photoassociation state is man- 
aged with a large intermediate detuning. To enable 
coherent phenomena, the two-photon Rabi frequency, 
%f2/<5, must be of the order of the two-photon spon- 
taneous decay lincwidth, il 2 j s /2S 2 , a condition that is 
not quite satisfied in the Wynar et al. [2£j experiments. 
On the other hand, atom-atom collisions are also an is- 
sue in achieving continuous-wave coherence, and avoid- 
ing them by increasing the laser intensity and decreas- 
ing the density could potentially bring rogue photodisso- 
ciation and laser-cloud overlap into play. Our opinion 
remains that coherent continuous-wave conversion can 
only be achieved by a difficult-to-impossible balance of 
laser intensity and condensate density. This bottleneck 
is probably not unique to the 87 Rb system. 

The present work has examined transient stimulated 
Raman adiabatic passage in photoassociation of a freely- 
interacting condensate in the presence of rogue and radia- 
tive decay. We have shown that rogue photodissociation 
to noncondensate modes is generally negligible, and that 
radiative decay, enhanced by collisional mean-field shifts, 
can be managed by adjusting the initial condensate den- 
sity and the atom-atom scattering length. In conclusion, 
a fully-optimized STIRAP scheme-i.e., where the detun- 
ings, pulse-sequence parameters (heights, widths, spac- 
ing, and number of sequences), particle number, trap 
parameters, and atomic scattering length are optimized- 
offers probably the best chance for achieving coherent 
atom-molecule conversion in photoassociation. 
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